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ABSTRACT 


Parameterizations of single nucleon emission from the electromagnetic interactions of cosmic 
rays with nuclei are presented. These parameterizations are based upon the most accurate 
theoretical calculations available today. When coupled with Strong interaction parameterizations, 
they should be very suitable for use in cosmic ray propagation through interstellar space, the 
Earth's atmosphere, lunar samples, meteorites and spacecraft walls. 



Introduction 


Galactic cosmic rays are very high energy particles confined to the region of our Milky Way 
galaxy. They consist of about 98% bare nuclei (stripped of all electrons) and about 2% electrons 
and protons (Simpson 1983). Of the nuclear component about 87% is hydrogen, about 12% is 
helium and the other 1% consists of heavier nuclei. Fe is the most abundant of these nuclei with a 
typical energy of about 1 GeV/N. Even though these heavy nuclei are not very abundant, they are 
very penetrating due to their large mass and high speed. 

An understanding of the interactions of galactic cosmic ray nuclei is important for several 
reasons: 

1) Knowledge of the cosmic ray spectrum at the top of the Earth's atmosphere and knowledge 
of the composition of the interstellar medium enables us to determine the cosmic ray spectrum at 
the source (Simpson 1983). 

2) Knowledge of the spectrum at the surface of the Earth and knowledge of the composition of 
the Earth's atmosphere enables us to determine the cosmic ray spectrum at the top of the 
atmosphere (Wilson, Townsend and Badavi 1987). 

3) The radiation environment inside a spacecraft, due to solar and galactic cosmic rays may be 
determined (Wilson and Townsend 1988). 

4) Studies of the history of extraterrestrial matter (such as lunar samples, meteorites and 
cosmic spherules and dust found in deep sea sediments) and also of the history of cosmic rays 
themselves can be made with knowledge of the production rate of various nuclides (Reedy 1987; 
Reedy, Arnold and Lai 1983). 

The basic nucleus-nucleus interaction that a cosmic ray undergoes can occur mainly via the 
Strong or Electromagnetic force. Strong interaction processes (Gyulassy 1981) have been studied 
extensively and quite recently the study of Electromagnetic processes in high energy nuclear 
collisions has begun (Bertulani and Baur 1988). 

In order to study the propagation of cosmic rays through interstellar space, the Earth's 
atmosphere or a spacecraft wall it is not enough to have only a good understanding of the nucleus- 



nucleus interaction mechanism. One must have an accurate theory of transport as well. Generally 
one uses a nucleus-nucleus interaction cross section as input to a transport computer code. These 
codes however can be veiy complex and therefore require simple expressions for the cross sections 
rather than the use of data bases or complicated theoretical expressions (Wilson and Townsend 
1988). Thus there has been a considerable effort to parameterize the cross section expressions so 
that the only required inputs are the nuclear energies and charge and mass numbers (Letaw, 
Silberberg and Tsao 1983; Silberberg and Tsao 1973; Townsend and Wilson 1986; Norbury, 
Cucinotta, Townsend and Wilson 1988; Wilson, Townsend and Badavi 1987). 

One approach to the parameterization of cross sections is to simply take all the available 
experimental data and fit a curve through it (Letaw, Silberberg and Tsao 1983; Silberberg and Tsao 
1973). Such an approach has certainly been useful and successful, but a much more satisfying 
parameterization would be one tied more directly to theory. It is the aim of the present work to 
obtain such a parameterization for the Electromagnetic (EM) part of the nucleus-nucleus interaction. 
One can then couple this with a similar theoretical parameterization of the Strong interaction 
process (Wilson, Townsend and Badavi 1987) to obtain a complete theoretical parameterization of 
the complete cross section. 

A preliminary parameterization of the EM process has already been presented (Norbury, 
Cucinotta, Townsend and Badavi 1988), which utilizes the Weizsacker-Williams (WW) method of 
virtual quanta (Bertulani and Baur 1988; Jackson 1975). However, since then the theory has been 
improved to include the effects of both electric dipole (El) and electric quadrupole (E2) interactions 
(Bertulani and Baur 1988; Norbury 1989a), which will henceforth be referred to as multipole 
theory in contrast to WW theory. These El and E2 effects modify the parameterization 
considerably. Also in the present work several different parameterizations are presented differing 
in degree of complexity. In addition much more data has become available with which to compare 
the parameterizations (Heckman and Lindstrom 1976; Olson, Berman, Greiner, Heckman, 
Lindstrom, Westfall and Crawford 1981; Mercier, Hill, Wohn, McCullough, Nieland, Winger, 
Howard, Renwick, Matheis and Smith 1986; Hill, Wohn, Winger and Smith 1988; Smith, Hill, 



Winger and Karol 1988; Hill, Wohn, Winger, Khayat, Leininger and Smith 1988; Hill, Wohn, 
Winger, Khayat, Mercier and Smith 1989; Norbury 1989b; Hill and Wohn 1989). The 
parameterizations to be presented below can then be combined with Strong interaction 
parameterizations such as the excellent parameterization by Wilson, Townsend and Badavi (1987). 
This combination should provide for much more accurate models of cosmic ray propagation 
through interstellar space, the Earth's atmosphere and spacecraft walls. 

The present work will only consider single nucleon emission from cosmic ray nuclei. This has 
been shown to be the dominant electromagnetic process. Other particle emission processes such as 
two-neutron emission have much smaller probability (Hill, Wohn, Winger, Khayat, Mercier and 
Smith 1989), and will be studied in future work. 


ELECTROMAGNETIC THEORY 


The EM theory has already been discussed extensively (Bertulani and Baur 1988; Norbury 
1989a) and only a few relevant details will be given here. The total nucleus-nucleus EM cross 
section is written as 
a = a E i + Oe 2 

= J[N E i (E) ofei (E) + N E2 (E) c E 2 (E)] dE (1) 

where N E i (E) is the virtual photon spectrum (of energy E) of a particular multipolarity due to the 
projectile nucleus and o E i (E) + a E2 (E) is the photonuclear reaction cross section of the target 
nucleus. (In principle the above equation should include other EM multipoles, but their effect is 
much less important.) A less exact expression is given by WW theory as 

CTww (E) = jNww (E) [g E i (E) + a E2 (E)] dE (2) 

where Nww (E) is the WW virtual photon spectrum. Bertulani and Baur (1988) have shown 
that 

Nww (E) = N E1 (E) = 1 1 Z 2 a -i- Ko - 1 £ (Kj 2 - Kq 2 )] (3a) 

bit p2 L 

N E2 (E)=l 1 Z 2 a^ [2 (1 - (3 2 ) Ki 2 + £ (2 - p 2 ) 2 KoKj-1 p 4 (Ki 2 - K 0 2 )] (3b) 


with 

^ _ Ebmin 
yp(1ic) 

where all of the Bessel functions K are functions of £. In the above equation E is the virtual 
photon energy, Z is the nuclear charge, a is the EM fine structure content, and b m ;n is the 
minimum impact parameter, below which the collision occurs via the Strong interaction. Also 



where c is the speed of light and v is the speed of the cosmic ray. 


P = ^- and 



The minimum impact parameter is given by 


bmin = Ro-i (T) + Ro-i (P) - d (5a) 

where Ro-i are the 10 per cent charge density radii of the projectile and target and d is an adjustable 
overlap parameter. An excellent approximation to Ro-i is (Norbury, Cucinotta, Townsend and 
Badavi 1988). 

Roi = (1.18 A 1/3 + 0.75) fm (5b) 

where A is the nuclear mass number. 

Jackson has provided high and low virtual photon energy approximations as 

NwwfEJ-llz^itln^-lp 2 ] (6a) 

for small and 


Nww(E) = lz 2 a-i-(l-i p 2 ) exp (- 2$) 
E p 2 2 


(6b) 


for large 

In equation (1) the El photonuclear cross section can be written in terms of the electric giant 
dipole resonance (GDR) cross section as 

o e 1 ( E ) = y — - 

l + [(E 2 -Ej DR r/E 2 lGDR ] 

where Egdr is the energy of the peak in the GDR cross section, Tgdr is the width of GDR, and 


rr - °TRK 
u m — 

tJgdr/2 

with the Thomas-Reiche-Kuhn cross section (Levinger 1960) given by 
Gt RK = 60 KZ MeV mb 

i\ 

where N and A are the neutron and mass numbers. The GDR energy is given by (Westfall, 
Wilson, Lindstrom, Crawford, Greiner and Heckman 1979) 

W2r2 


( 8 ) 


(9) 


Egdr =*c [ - c R ° (1 + u - 1 + g - + 3u -e)] 
8J 1 + e + u 


- 1/2 


GO) 



with 


u=21a-W (11) 

and 

Ro =r 0 A 1/3 (12) 

where e = 0.0768, Q' = 17 MeV, J = 36.8 MeV, ro = 1.18 fm, and m* is 7/10 of the nucleon 
mass. Note that other expressions for Egdr such as 80A _1/3 (Bertulani and Baur 1988) provide 
very inaccurate results for light nuclei. Equation (10) is accurate for all mass regions. 

The E2 cross section is dominated by the giant quadrupole resonance (GQR). The main 
contribution to single nucleon emission (Bertulani and Baur 1988) comes from the isoscalar 
component given by (Bertrand 1976) 


o E2 (E) = OewsrE 2 

1 + (E 2 - Eqqr) 2 / E 2 !' 2 GQR 

with the energy- weighted sum rule (EWSR) cross section 

0.22 ZA 2 ^ [Lb MeV 4 

CtEWSR = 1 

n Tgqr/2 

where f is the fractional exhaustion of the EWSR (Bertrand 1976) and 


(13a) 


(13b) 


Egqr = ^|- (14) 

Finally, all of the above cross sections refer to total absorption cross sections. To obtain the 
reaction for proton or neutron emission they must be multiplied by the proton or neutron branching 
ratios. The proton branching ratio has been parameterized by Westfall et al (Westfall, Wilson, 
Lindstrom, Crawford, Greiner and Heckman 1979) as 

g p = Min [Z/A, 1.95 exp (-0.075 Z)] (J 5 a ) 

where Z is the number of protons and the minimum value of the two quantities in square brackets 
is to be taken. Assuming that only single nucleon emission occurs, the neutron branching ratio is 
gn = 1 - g P (15b) 



Comparison Between Theory and Experiment 


In the above paragraphs I have provided the basic equations to be used in the present work. 
However in analyzing the validity of the basic EM theory one uses only equations (3) - (5) and 
instead of equations (7) - (15) one uses actual experimental data for the photonuclear cross 
sections. A detailed study of the validity of this EM theory has been made (Norbury 1989a, b, c, 
d) and the results from this work are presented in Table 1, both for the WW theory and the separate 
El and E2 multipole theory calculations, and are compared to experimental data. A detailed 
discussion is to be found in Norbury 1989a, b, c, d, but the following features are to be noticed. 
Both WW and multipole theory give reasonably good results although multipole theory is 
somewhat better. It is found that electric quadrupole (E2) effects are not significant for proton and 
neutron emission from 12 C, 16 0 or 18 0. However, E2 contributions are substantial for neutron 
emission from 59 Co, 89 Y and 197 Au, generally leading to improved agreement between theory and 
experiment Notable disagreements occur for 139 La projectiles (1.26 GeV/N) where the theoretical 
Oei + Oe 2 are too big. Quadrupole effects improve the theoretical results for 16 0 projectiles at 60 
and 200 GeV/N, although the theoretical cross sections are still too small. In general it has been 
found (Norbuiy 1989a, d) that electric quadrupole effects are an important component in nucleus- 
nucleus collisions and that these effects can be calculated accurately. 



Parameterizations 


As mentioned above in testing the basic WW and multipole EM theory one uses experimental 
data for the photonuclear cross sections. However this is not a practical procedure for use in 
cosmic ray transport codes and instead my approach will be to use expressions (7) - (15). 

In the present work I shall discuss three separate parameterizations of the above EM theory for 
use in cosmic ray transport codes. These will be presented in decreasing order of accuracy, but the 
aim is to provide parameterizations that will be useful in different contexts. 

Parameterization #1 of Multipole Theory 
This is the most accurate parameterization and uses the following equations: 

1) Equation (1) is used for the total nucleus-nucleus EM cross section. The integration is 
done numerically using the Trapezoidal Rule. 

2) Equations (3) are used for the virtual photon spectra Nei (E) and Ne 2 (E). 

3) Equations (5) are used for the minimum impact parameter with the overlap parameter 
adjusted to give the best fit to data at d = -1.5 fm. 

4) Equations (7) - (14) are used for the photonuclear cross sections. 

5) The width Tgdr in equations (7) and (8) is set at 

T G dr= lOMeV for A <50 

= 4.5 MeV for A > 50 (16) 

and Tgqr in equations (13) is set at 

Tgqr = 2.5 MeV for A > 180 

= 4.5 MeV for 70<A<180 
= 5.5 MeV for 19 < A <70 

= 3.0 MeV for A < 19 (17) 

These values for Tgrd are discussed in Norbury, Cucinotta, Townsend and Badavi 
(1988) and for Tgqr in Bertrand (1976). 

6) The fractional exhaustion of the Energy-Weighted Sum Rule in equation (13b) is given 
by (Bertrand 1976) 

f = 0.9 for A> 100 
= 0.6 for 40 < A < 100 
= 0.3 for 40 < A 

7) The proton and neutron branching ratios are given by equations (15). 


( 18 ) 



The results of the above parameterizations are given in Table 1. It can be seen that it agrees 
extremely well with the multipole theory. Thus I regard this parameterization #1 of the multipole 
theory as describing very accurately the most advanced state-of-the-art EM theory. Agreement 
between this parameterization and experiment is, of course, of the same quality as between the 
multipole theory and experiment. 

Parameterization #1 of WW Theory 

WW theory gives a simpler treatment of the virtual photon field and is included here for the 
sake of completeness. The only difference between parameterization #1 of WW theory and 
parameterization #1 of multipole theory is that equation (2) is used for the total cross section 
instead of equation (1). Results are listed in Table 1 and are fairly comparable to the 
parameterization #1 of the multipole theory. 

Parameterization #2 of Multipole Theory 

A difficulty that might occur in some cosmic ray transport theories is the necessity of having to 
do a numerical integration in equation (1) every time a is to be evaluated. To get around this, 
parameterization #2 is based on the technique of Bertulani and Baur (1988). This involves taking 
Nfii (E) outside of the integral in equation (1) and evaluating Ngi (E) at Eg dr (see equation 10) and 
Ne 2 (E) at Egqr (equation 14). The remaining integral is evaluated from sum rules. That is 
(Bertulani and Baur 1988), equation (1) becomes 

<J - Nei (Egdr) J oti (E) dE + N E2 (Egqr) eJq R J^^dE (19) 

with the sum rules 

j Obi (E) dE = mb MeV ( 20a) 

A 

and 

J °E2_® dE =f QJSJZA PL mb MeV- 1 (20b) 

g2 1000 

Bertulani and Baur (1988) claim that this is an accurate procedure. However, I found it necessary 
to change d to d = -2.4 fm (see equation 5a) in order to give good comparison to experiment. 


In the present parameterization #2 of multipole theory items 1) - 3) of parameterization #1 were 
changed to those discussed in the preceding paragraph. Note especially that a numerical integration 
is no longer necessary. Items 4) - 5) are no longer relevant. Items 6) - 7) remained the same. 
Results are again listed in Table 1. With the new value of d = -2.4 fm parameterization #2 agrees 
well with parameterization #1 (which used d = -1.5 fm). 

Parameterization #2 of WW Theory 

WW theory is again included for completeness. In this case equation (2) was replaced with 

Cww ~ Nww (Egdr) J o E i (E) dE + Nww (Egqr) Egq R dE (21) 

with the same sum rules in equations (20). Results are listed in Table 1. 

Parameterization #3 

Parameterizations #1 and #2 require the evaluation of Bessel functions as indicated in equations 
(3) for the virtual photon spectra. In the interest of providing an even simpler parameterization that 
could be used on a pocket calculator for rough estimates of the cross section, a third 
parameterization is presented. The E2 cross section was ignored and equations (1) or (2) were 
replaced with 

a = Nww (Egdr) J ofei (E) dE (22) 

Note that this is identical to neglecting the GQR in equation (1). The sum rule in equation (20a) 
was used for the integral. Nww (Egdr) was evaluated using equations (6), with (6a) used for 
£, < 0.5 and (6b) for £, > 0.5. This prescription avoids the evaluation of Bessel functions and 
almost allows one to calculate a in one's head. In this case the value of d was d = +1.0 fm. Items 
4) - 6) are not relevant and item 7) was again used. The results are presented in Table 1 and are 
seen to give surprisingly similar results to the other parameterizations. 




As discussed in previous work (Norbury 1989a, d) the multipole theory is generally more 
accurate than WW theory. This is also true for the above parametrizations as can be seen from 
Table 1. 

However, WW theory and multipole theory do not describe 18 0 very well, and the 
parameterizations are even worse. I trace this to the fact that the branching ratio equations (15) do 
not work well for nuclei off the stability curve. 

Both WW theory and multipole theory do not describe 197 Au very well either, but the 
parameterizations do a somewhat better job due to the choice of the overlap parameter d. There 
seems to be a problem also for veiy high energies especially 200 GeV/N. 

Apart from these problems the multipole theory and multipole parameterizations (#1, #2 and 
#3) seem to describe the data quite accurately. 

As regards which parameterization to use, they all seem to do an equivalent job in describing 
the data. This of course is because a different value for d was chosen for each. Even the 
parameterization #3 does quite well, although it is a little high for nucleon emission from the lighter 
nuclei. 

Given the above problems with 18 0, 197 Au and 200 GeV/N I recommend that the above 
parameterizations be used 0 only with nuclei on the stability curve, “) for nuclei lighter than 197 Au 
and *“) for energies less than 10 GeV/N. These requirements should not be too restrictive in 
Cosmic Ray work because most nuclei have energies of around 1 GeV/N and the most abundant 
nuclei are not much heavier that 56 Fe (Simpson 1983). Having to deal only with nuclei on the 
stability curve is probably the most severe restriction. 

Parameterizations #1, #2, #3 decrease in order of accuracy, but, as discussed above, not by 
very much. I would recommend using the most accurate parameterization (#1), but if one's 
computer codes are such that it would save CPU time by using either #2 or #3, then I would 
recommend their use. However, one should perhaps be careful about using parameterization #3 


for light nuclei. I recommend the multipole parameterizations, but I do not recommend the use of 
the WW parameterizations. 

Finally, by combining the above EM parameterizations with the Strong Interaction 
parameterization of Wilson, Townsend and Badavi (1987), which is not subject to the same 
restrictions as above, transport of cosmic rays through matter can be described very accurately. 
Future work will involve parameterization of both multiple nucleon emission (a much smaller 
effect) and also neutron branching ratios for nuclei off the stability curve. 
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